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Abstract 

Starting from critical RSOS lattice models with appropriate inhomogeneities, 
we derive two component nonlinear integral equations to describe the finite 
volume ground state energy of the massive (p^ ^ a( jj perturbation of the 
SU(2)k x SU(2)y I 'S77(2)fc + fc/ coset models. When k' — > oo while the value 
of k is fixed, the equations correspond to the current-current perturbation 
of the SU(2)k WZW model. Then modifying one of the kernel functions of 
these equations, we propose two component nonlinear integral equations for 
the fractional supersymmetric sine-Gordon models. The lattice versions of our 
equations describe the finite size effects in the corresponding lattice models, 
namely in the critical RSOS(/c,(/) models, in the isotropic higher-spin vertex 
models, and in the anisotropic higher-spin vertex models. Numerical and an- 
alytical checks are also performed to confirm the correctness of our equations. 
These type of equations make it easier to treat the excited state problem. 



1 Introduction 



The massive id ac |j perturbation of the SU(2) k x SU(2) k > / SU(2) k+k > models 
belong to an interesting set of integrable quantum field theories (QFT). In the 
following we will denote them by -MA^ k+k , +2 - This family of models includes as 
particular cases some other interesting QFTs like the asymptotically free current- 
current perturbation of the SU(2) k WZW model and the 0(4) nonlinear sigma 
(NLS) model. The spectrum and interactions are governed by two dual fractional 
symmetries and Q^ k \ one of which turns into the fractional supersymmetry 
for k = 1 pi E] • The spectrum of the theory consists of a set of massive relativistic 
kinks of the same mass P] Their factorised scattering was conjectured in pQ. 
The kinks K aa i^y separate different degenerate vacua, each vacuum being labeled 
by two indices aa! . a and a' run independently over the states of nodes of an A k+ \ 
and A k i + i Dynkin diagrams respectively, a e {0,1,..., A;} , a' G {0, 1, &;'} . The 
stable kinks K aa i y w follow the Dynkin diagram adjacency rules a — 6± 1, a' — 6'± 1 
and the kink-kink scattering matrix appears as the tensor product 

s{e) = sf SG {e)®s^ G {e), (i) 

where Sr SG (9) is the scattering matrix of the massive 0i3 perturbation of the min- 
imal model SU {2)i x SU(2) k / SU(2)i +k , which is an integrable restriction of the 
sine- Gordon model 

The fractional supersymmetric sine-Gordon models (FSSG) were first proposed 
as the models obtained by "unrestricting" the S-matrices ((TJ) describing M.A~l k+k , +2 
j2j. Thus the scattering matrix of the FSSG models is of the form 

S(0) = S ( *l G (6)®S^ G (0), (2) 

where S^ G (6) is the S-matrix of the sine-Gordon model at coupling (5 The 
particle content of the FSSG consist of a soliton and an antisoliton of mass m 
together with a number of breathers of mass 

M 3 = 2m sin f^f] , j = 1, ... < 1/f (3) 



where 

S l/fc-/3 2 /87r' 1 } 
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and there are no breathers in the repulsive regime £ < 1. 

Once the spectrum and the associated S-matrix is known, the TBA method can 
be employed for calculating the ground state energy of the model in finite volume. 
The TBA equations of the A4A^ r model was conjectured in 6J and the equations 
can be encoded in an A r _3 Dynkin-diagram (see figure^.) The TBA equations 
of the FSSG models were proposed in [7], and at the special ^ = ^-j-, 1 < r G N 
values of the coupling, it can be encoded into a T> T Dynkin-diagram (see figure 
HJ>), but for general values of the coupling, the TBA integral equations take much 
complicated form, and the number of unknown functions of the equations depend 
on the continued fraction form of r The solutions of the TBA equations also 
satisfy the well-known Y-system equations jH]. 

r 

O O O • O O O a o — c^d o — • — c o — / b 

1 2 k-1 k k+1 r-4 r-3 1 2 ^ k k+1 r \_ 

r-l 

Figure 1: Graphical representation of Y-system of the a.) A4A^ r and b.) 
FSSG model 

In |Zj it was noted that if a model has an S-matrix in the form of a direct product 
Sg ® Sh, and the TBA equations are already known for the models described by 
S-matrices So and Sh, and the individual TBA equations are encoded on Dynkin- 
like diagrams of type G and H respectively, each having one massive node, then the 
TBA equations for the model with the direct product S-matrix can be obtained by 
gluing the individual TBA equations together at the massive node. This method 
can also be applied for our perturbed coset model due to the tensor product form 
of the S-matrix. 

Another approach for calculating finite size effects in a QFT is the use of some 
integrable lattice regularization of the QFT which can be solved by the Bethe Ansatz 
method. In models that can be solved by Bethe Ansatz, the well-known T- and Y- 
systems [TIH ITT] naturally appear and from them, using the analytical properties 
of the T-functions, one can easily derive the same TBA equations which can be 
obtained from the S-matrix. Using such an integrable regularization, an alternative 
method was worked out by Destri and de Vega ^2] for calculating finite size effects 
in the SG model. This is called the nonlinear integral equation (NLIE) technique. 
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The advantage of this method is that one gets a single integral equation for any real 
value of the coupling constant, and the method can be extended to excited states, 
too El E] • The problem is that this method can only be applied when the 
ground state of the model is formed by real Bethe Ansatz roots. In most cases this 
is not true. For example the ground state of the spin-S* XXX-chain is formed by 
2S-strings PH El El E] • 

Recently J. Suzuki [20J managed to derive some new types of nonlinear integral 
equations for describing the thermodynamics of the higher spin XXX-model at 
finite temperature, which can be regarded as the generalization of the spin-^ case 
which was treated successfully earlier [121 EI] • These new equations can be regarded 
as a particular mixture of TBA and NLIE, the special case of which is the Destri-de 
Vega equation [TJ] in the spin-^ case. 

Motivated by the gluing idea [Zj and Suzuki's results [20] Dunning proposed 
nonlinear integral equations [22] similar to Suzuki's equations for describing the 
finite size effects in some perturbed conformal field theories whose S-matrices can be 
written in the form of a direct product, comprising the AiA~^ r models, the fractional 
super symmetric sine-Gordon models (FSSG), and the SS-model respectively. 

These equations can also be regarded as a mixture of TBA and NLIE, but these 
equations have the advantage that the number of unknown functions is less than 
it is in the TBA equations, hence these equations are much more convenient for 
numerical studies. Another advantage is that if the model has a single coupling 
constant then the equations are valid for all real value of the coupling constant. 

Later two component nonlinear integral equations were derived for the ground 
state energy of the 0(4) nonlinear a-model [2S] in the context of the light-cone 
approach [2UI22], formulating the 0(4) nonlinear cr-model as the S — > oo limit of 
the inhomogeneous spin-S vertex model, and by modifying the kernel functions of 
these equations similar nonlinear integral equations were proposed for the finite size 
effects in the SS-model for finite values of the couplings |2"3~j . 

In this paper, formulating MA^ r model as the appropriate continuum limit 
of the critical, inhomogeneous RSOS(fc,g) lattice model, we derive two component 
nonlinear integral equations for describing the ground state energy of the model. The 
lattice version of these equations describe the finite size effects in the corresponding 
critical RSOS(fc,g) lattice model. 
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When the allowed values of the heights tends to infinity (r — ► oo), we obtain the 
equations corresponding to the current- current perturbation of the SU(2)k WZW 
model. The lattice version of these equations describe the finite size effects in the 
isotropic spin-/c/2 vertex model. 

Then formally releasing the restriction in the two component nonlinear integral 
equations of the M.A\ r model, we propose two component equations for the FSSG 
models, the lattice version of which correspond to the case of the anisotropic higher 
spin vertex model. 

Due to the minimal number of components, our new nonlinear integral equations 
are more convenient for numerical studies than the corresponding TBA equations. 
Moreover, the minimality of the number of components makes the treatment of the 
excited state problem much easier than in the context of the TBA approach. 

The paper is organized as follows. In section 2, we recall the main results of the 
light-cone lattice approach to the M.A\ r models. In section 3, we summarize the 
most important properties of the critical RSOS(A;,g) lattice models, comprising the 
T- and Y-systems, which follow from the Bethe Ansatz solution of the model. In 
section 4, using the light-cone lattice approach and Suzuki's method we derive 
nonlinear integral equations for describing the finite size effects. In section 5, for- 
mulating the RSOS(/c, q) lattice model as an RSOS(A;', q') model with k' = r — k — 2, 
q' = r — q — 2, we derive another nonlinear integral equations for the same model. In 
section 6, we glue together the two nonlinear integral equations to a much simpler 
equation. This equation describes the finite-volume ground state energy of our per- 
turbed coset model. In section 7, we perform some analytical and numerical tests on 
our two component equations. In section 8, we propose nonlinear integral equations 
for describing the finite- volume ground state energy of the fractional supersymmetric 
sine-Gordon models by formally releasing the restriction in the equations obtained 
for the A4A^ r models. In section 9, we briefly summarize how the two component 
nonlinear integral equations look like in the case of the solvable lattice models. The 
summary and conclusions of this paper are given in section 10. 



4 



2 The light-cone lattice approach to the MA£ r 
models 

In the following we use the integrable light-cone lattice regularization of our model 
j2l)] I24j . Starting from the inhomogeneous critical RSOS lattice model, in a 
local integrable lattice regularization of the Hamiltonian of the AiA~^ r model was 
proposed. As a result of the locality of the lattice Hamiltonian, in the continuum 
limit of this regularization, one obtains the A4A^ r model, plus a decoupled SU (2) k x 
SU X2) r _fc_ 2 / 'SU (2) r _ 2 coset conformal field theory. In the following to avoid such 
"particle doublings", following the lines of [21], we define a nonlocal Hamiltonian 
for the A4A^ r model from the same lattice model. The final result of the two 
regularizations are the same in the continuum. The fields of the regularized theory 
are defined at sites ( "events" ) of a light-cone lattice and the dynamics of the system 
is defined by translations in the left and right light-cone directions. These are given 
by inhomogeneous transfer matrices of the RSOS(/c,g) model considered on the 
critical line separating the regimes III and IV. 

This approach is particularly useful for calculating the finite size dependence of 
physical quantities. We take N points (N is even) in the spatial direction and use 
periodic boundary conditions. The lattice spacing is related to I, the (dimensionful) 
size of the system : 

a = AT ( 5 ) 
The energy (E) and momentum (P) of the physical states can be obtained from 
the eigenvalues of the light-cone transfer matrices (j!4jl : 



j&(h+p) _ ^fc ) ( x o + %%[H-P) _ \~ x o ~ i) 

~ T (fc) (x + *(£: + l))' " T (fc) (-^o-*(£: + l)) 



Besides the usual procedure, taking the thermodynamic limit (N — > oo) first, fol- 
lowed by the continuum limit (a — > 0) one can also study continuum limit in finite 
volume by taking N — > oo and tuning the inhomogeneity parameter x$ simultane- 
ously as 

2 , 4 2 , 2N 
x = - log = - log — , (7) 

7r ma 7T ml 

where the mass parameter m is the infinite volume mass gap of the theory. If we take 
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this continuum limit we get the energy and momentum eigenvalues of the M.A~l 
models. 



3 RSOS(fc,g) lattice models 



The restricted solid-on-solid models of Andrews, Baxter and Forrester j2Zj are FRF 
or interaction-round-a-face models [2B]- The heights or spins a,b,c,d, etc. at each 
site of the square lattice take the values l,2,...,r-l. These spins are subject to the 
nearest neighbour constraint a — b = ±1 for each pair of adjacent spins a, b. The 
spins therefore take values on the Dynkin diagram of the classical Lie algebra v4 r _i. 
The statistical weight assigned to an elementary face of the lattice is zero unless 
the spins around the face satisfy the adjacency rules described later. The weights 
of allowed faces are given by 
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sin(A) 
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S a = sin(aA). 



(_l)(a-c)/2 (8) 



(9) 



The spectral parameter u is related to the spatial anisotropy, A is the crossing 
parameter and 5 is the Kronecker delta. 

The critical RSOS(A;,g) lattice models are obtained by fusing k x q blocks of 
face weights together. In the following we consider the model on the regime III/IV 
critical line. The face weights of the fused RSOS(fc, q) models, with max(/c, q) < r, 
are obtained by the fusion process |2H1- The RSOS(/;;,g) face weights are given 
explicitly by [23 E01 



W k ' q 



u 



a q+l b q+ \ 

ai bi 



q-2 q-1 

nn 



sin(A) 



n=0 j=0 
a n+l b n+ i 

Q>n. bn 



sin(w + (n — j)X) 
u + (k- 1)A 



(10) 



ai,...,a q n=l 

independent of the values of the edge spins b 2 , .., b q , where the fcxl face weights are 
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given in turn by 



W 



k.l 





bk+i 
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(22,...,flfe 71=1 



On O^+l 



u + (n — k)X 
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independent of the values of the edge spins &2>--,0fc- For the RSOS(A;,g) models, 
adjacent spins or heights are subject to constraints 

< (oj — aj + m)/2 < m, m < + a.,- < 2r — m, (12) 

where m is equal to k for a horizontal pair and g for a vertical pair. 

The RSOS(fc, q) lattice models are exactly solvable. In particular, the fused face 
weights satisfy the generalized Yang-Baxter equation 




u + v \ W 



u + v\ W 




(13) 



This is an immediate consequence of the elementary Yang-Baxter equation satisfied 
by the lxl face weights and leads to commuting transfer matrices. Suppose that a 
and b are allowed spin configurations of two consecutive rows of an iV column lattice 
with periodic boundary conditions. Assuming that N is even the matrix elements 
of the RSOS(/c,g) transfer matrix with alternating inhomogeneities are given by 

N 



i\TW(x)\b) = Y[\V k ' q 



bj b j+1 

Oj CLj+l 



ITT 



(x - (-l) j+1 x ) k + 1 



Q 



2r 



A 



(14) 



where a N+ i = a x and b N+1 = 6 X , and x is the inhomogeneity, which is zero if one is 
interested in the finite size dependence of the energy levels of the lattice model, and 
must be tuned according to (JZJ) to obtain a continuum field theory in finite volume. 
For a fixed value of k the Yang-Baxter equations imply the commutation relations 



Ti k \x)T^(y)=T^ (y)T$ k \x). 



(15) 



The eigenvalues of these mutually commuting transfer matrices satisfy the so-called 
T-system functional relations [3T] 

T} h \x - i) T«(x + i) = /f (x) + T^(x) T^{x\ (16) 
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where 

/« = T (fe) (x + i(p + 1)) T (fc) (* - i(p + 1)), (17) 

and 

k-l 

T {k } = 0, T (fc) (x) = H<j>(x + i(k - 1 - 2j)), (18) 

i=o 

sinh^d (z - x )) sinh^d (x - x )) 
sin N (X) 

These T-system equations are the same as the ones of the anisotropic higher-spin 
vertex model [3U]- The RSOS(A;,g) transfer matrices satisfy the symmetry [3U] 



4>(x) = ° y ^ 7vr; 2 — . (19) 



Tf \x) = YTi%(x -ir), q = 1, 2, r - 3, (20) 



r-2-q 

where Y is the height reflection operator 



N 



(a|Y|b) = nC J , [Tf (*),Y] = 0. (21) 

3=1 

It follows that the fusion hierarchy closes with T r ^(x) = 0. The functions Tq(x) 
are all periodic on the complex plane with period 2ir 

T^ k) {x)=T^ k \x + 2ir). (22) 

It will be important for our later considerations that it can be shown [HD! that 
the transfer matrices Tq k \x) (pPfj) of an RSOS(A;,g) model differs from the transfer 
matrices T\v\x) of an RSOS(&', q') model only in some trivial normalizations, if 
k' = r — k — 2, and q' = r — q — 2: 

Tf(i) ~Tp{x), k' = r-k-2, q' = r-2-q. (23) 

The solution of the T-system (jlfij) can be characterized by a set of Bethe roots 
{xj,j = 1, kN/2 }, which are the solutions of the Bethe Ansatz equations 

T?\x j + i) ^ QW(x j + 2i) 



where 

kN/2 



1 / \ \ 

Q {k \x)= JJsinh ^-(x- Xj )j 



(25) 



e luJT = -(-l) kN/2 Y, u^O. (26) 

where Y — ±1 is the eigenvalue of the operator Y. The eigenvalues of the transfer 
matrices (fT4"j) are of the form [3U1 EH] 

p+i 

T^) = £A! p) (x), (27) 
i=i 

where 

A, {X)-l [X + l{2l p I)) Q{k){x + t{2l _ p _ 1))Q{k){x + t{2l _ p _ 3) y 

(28) 

From a T-system (jlfij) one can define a Y-system as follows 



(k)r \ ^j-l( x )^j+l( x ) ,, )fV . 



Y} , (x) = l + y) >(x) = -1 -^f- . (30) 

fj ( x ) 

These functions satisfy the Y-system equations EDI EE] 

yf\x + i)yf\x - i) = Y^Y^x), j = 1, ...r - 3. (31) 

It follows from (J2Uj) that the Y-system (|29H31|) also closes (ie. yfH 2 {x) = 0). One 
can redefine the T-system elements with 

T j k \ x )^f^(x) = a J (x)Tj k \x), (32) 

where (Tj(x) satisfies the relation 

<jj(x + i)<jj(x — i) = aj-i(x)o-j + i(x), (33) 

then the new (x) functions satisfy the same T-system relations as Qlfij) but with 
different ff\x) functions: 

f?\ x ) - ff\x) = a 3 {x + i)a {x - i)ff\x). (34) 



9 



This is called gauge transformation. Under such a transformation the Y-system 
(I29H31)) is invariant. As we will see later the energy of the model can be expressed 
by an element of the gauge invariant Y-system (see eq. ()104|) later). On the other 
hand one has the freedom of choosing that gauge for the T-system (j!6|) which is 
convenient for the particular calculations. For later use it is important to note that 
the Y-system f|29H3 lj) of an RSOS(fc, q) model following from (|23~j) is identical to the 
one of an RSOS(fc', q') model ie. 

y( k \x) = y { f\x), k' = r-k-2 q' = r-q-2. (35) 

From this point of view an RSOS(fc, q) model is equivalent to an RSOS(A; / , q') model, 
and so in any particular calculation we can choose the formulation which is more 
convenient. 



4 Nonlinear integral equations I. 

In this section we apply Suzuki's method |2H| to obtain a set of nonlinear integral 
equations for the ground state energy of our model in finite volume. Consider the 
following gauge transformation of the T-system (|T6*j) 

f f\x + i)T}*>(x -i) = fj\x) + fl k \{x)fM x (x), (36) 

where 

p 

fi k \x) = H<P(x + t(p-k- 2j)) <p(x -i(p-k- 2j)), (37) 
i=i 

and 

fl fc 1 ) (.x) = 0, f (fc) (x) = l. (38) 

The solutions of these equations are of the form j2[ 



P+i 
i=i 



where 



x(p) M = ^(p+i-v) 7<p), x qw(s + »(p + i)) gw(x-»(p + i)) 

' { ) Wl { ) Q( k )(x + i(2l-p-l)) QW(x + i(2l-p-3)Y 



(40) 
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p-l+1 l-l 
^\ p \x) = ] [ 4>(x + i(p-k- 2j + 1)) Y[4>(x -i(p-k- 2j + 1)). (41) 

3=1 3=1 

We define the following auxiliary functions 



yfw = S|giM, > = 1 * (42) 

if'H-'^fH- ^lr 11 . i = i.-.* («) 



Ki) = Ai»(x + i) + - + y»(x +< ) | 6(x) = l + (44) 



tt»)= ^ (j "Uy + ?' (j " <) - «<*> = !+«*>■ («) 
Aj (x — z) 

From (|36H4H) it follows that the auxiliary functions ()42M5|) satisfy the following 
functional relations: 

ft% + z) = e-*» J] 0(x + 2y) |^|| B(x), (46) 

^ (fc) (- - = ^ fl «K* - «i) %rg ( 47 ) 

- b(x) _ e -*<* + u« <M g (fc) (g-»fc-2Q f (fe) f v , . 

= F fc W (a;), (50) 

yf\x + z)yf } (z - z) = Y}1\(x)Y}1\(x) j = 1, . . . , fc - 2, (51) 

v&C* + Oi/ffiO* + «•) = (52) 



n 



fi^x + i)f^(x -i) = f { k k \(x)Y£\(x). (53) 
We introduce two other auxiliary functions 

q[ k \x) = Q( k \x-tk), ^2 ] {x) = Q {k \x + ik) = ¥ k \x-2i{r- k)). (54) 

In order to be able to derive integral equations from these functional relations one 
needs to know the positions of the zeroes and the poles of the auxiliary functions 
(J42H45)) . Due to the relations (|4fifl5Hj) we only need to know the zeroes of 
and Tj (x). For the ground state of our model the zeroes of Q^ k \x) form y pieces 
of /c-strings, but for finite N values there are deviations from the string hypothesis 
|32| . From [22] one can see that Bethe roots having the largest deviations from the 
imaginary positions prescribed by the string hypothesis are those that would have 
imaginary parts ±(k — 1) according to the string hypothesis. These deviations are 
always less than 1/2 and this is important because the analytic properties of (x) 
(^^(x)) are influenced mainly by these roots on the upper (lower) half plane of the 
complex plane near the real axis. 

The transfer matrices fj(x) j = 1, . . . , k (I36|) have no zeroes in the ground 
state in the "main" strip < \lmx\ < 1. We can list the strips where the auxiliary 
functions are analytic and non zero (ANZ). 



y[ k \x) 


ANZ 


-2(r - k) - 1/2 < Im x < 1/2, 


mo 


d 2r 






¥ k \x) 


ANZ 


-1/2 < Im x < 2(r - k) + 1/2, 


mo 


d 2r 






b(x),B(x) 


ANZ 


< |Im x\ < 1/2, mod 2r 










h(x),B{x) 


ANZ 


< Im x\ < 1/2, mod 2r 








(55) 


yf\x) 


ANZ 


< Im x\ < 1, mod 2r j = 


1,. 


..,k- 


1, 




Yf k \x) 


ANZ 


< |Im x\ < e, mod 2r j = 


h- 


..,k- 


1 


e > 0, 


7f } (x) 


ANZ 


< Im x\ < 1, mod 2r j = 


1,- 


..,k. 







We introduce new variables by shifting the arguments of b(x),B(x) and b(x),B(x) 
by ±17 |2D] 



a (x) = b(x — ij), U (x) = B{x — ij) = 1 + a (x), (56) 
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ao(x) = b(x + ij), Uq{x) = B(x + ij) = 1 + a (x), (57) 

where 0<7< 1/2 is an arbitrary real parameter. This shift is necessary because the 
original functions b(x),b(x) have zeroes on the real axis. Due to the ANZ property 
of T^\x) (JH5j) and the fact that \im.\ x \^ 00 -^2\ogT < ^\x) = the following relation 
holds 



= J dx log if \x - i) e iqix - i] - J dx logTf \x + i) e iq{x+i \ (58) 

— oo — oo 

From ()46|) and ()47|1 one can express (x ± z) with the auxiliary functions and by 
substituting these expressions into (JH%|) one gets in Fourier space 



dHm\ k \q) = TrNqe { 



sinh(2g) sinh(rg) 



e (r-fc+l-7) 9 cPlU (q) e (r-k-l+y)g d 2 lU (q) 

_l_ ( o9 ) 

4 cosh(g) sinh[(r — k)q) 4 cosh(g) sinh[(r — k)q] ' 



dmf>{q) = nNqe 



sinh(2g) sinh(rg) 



+ e -(r-fc+l+ 7 )g d?lU Q (q) e -fr-fc-l-7)q g^(g) q 

4 cosh(g) sinh[(r — k)q] 4 cosh(g) sinh[(r — k)q] ' 



where we introduced the notation 



oo 



d~ 2 lF(q)= I dx e iqx -J-logF(z). (61) 
V era? 

— oo 

One can derive similar relations for d~ 2 li/j k \q)'s and cPIYj (?)'s from (|51)L and 
cPlT^^q) and cPlY^^q) from (jHHJl- Substituting these relations into the defini- 
tions of cto( a; ) an d Qo(^)) one obtains k + 1 algebraic relations in Fourier space. 
After taking the inverse Fourier transformation of these relations and integrating 
twice over x determining the integration constants from the x — > oo asymptotics 
of the auxiliary functions and exploiting that following from (J2fijl u = ir/r for the 
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ground state, we get Suzuki's "half TBA-NLIE" type of equations 



\ogy[ k \x) 
\ogyf\x) 
log 2/^(2) 
log a Q (x) 



(K *\ogY 2 (k) )(x 



(K*iogy ( 5) + (ir*iogy^)( 



(*)• 



X 



j = 2,...,k-2, 



(K * log y£ } 2 ) (x) + (K + t * log U )(x) + (A~ -* * log C/ ) (x) , 
D^x - 27) + (G * \ogU )(x) - (G +2 ^ * \ogU )(x) 



(62) 



+ (K-~<*\ogY£\)(x)+i 



7T 



r — k 



— i57T, 



loga (x) = V N (x + t 1 ) + (G*logU )(x)-(G- 2 ^*logU )(x) 



(fc)w . . 7T 



.T 



r — k 



where (K*f)(x) = J dy K(x — y) f(y) is the convolution and the "source" function 
T> N (x) on the lattice reads as: 



T>n(x) = ziVarctan 



sinh 



cosh (2fa) 



and 



5 = (iV/2) mod 2 G {0, 1} 



X ° = ^ l0g (^)' (63) 



(64) 



the value of which is zero (5 = 0) for the ground state of the AiA^ r model. The 
kernel functions of (152*1) are of the form 



K(x) 



1 



4 cosh(7ra;/2) 



(65) 



G{x) 



dq ^ iqx sinh[(r - k - l)g] 
2tt 2 cosh(g) sinh[(r — k)q] 



and we have used the notation 



f ± »{x) = f(x±iri). 



(66) 



(67) 



We get the continuum limit of eqs. (J"**2*|) by taking the N 00 limit. The ker- 
nel functions do not change because they are independent of N, but the "source" 
function changes in the continuum as 



V(x) 



lim Vn(x) 



-ml cosh ^— j . 



(6* 
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In the TBA language the complex auxiliary functions ciq(x) and clq(x) resum the 
contributions of those Y-system elements whose index is larger than k — 1. Equations 
(|62|) are graphically represented in figure Eh- 



RSOS(fc,<?) — ► \ / \ ' — RSOS ( fc '.9') 

o — o o — • — o O O J a ( o — o q — • — o o — o 

1 2 k-1 \k k+1 r-4 r-3/ V" 3 r ~ 4 k '+V k' k'-l 2 1 



a (x),a (x),G(x) a(x),a(x),G(x) 

Figure 2: Graphical representation of nonlinear integral equations a.) I and 
b.) II. 

The big "bubble" denotes the complex auxiliary functions which resum the con- 
tributions of those TBA nodes, which are inside it. In our notation the names of the 
complex unknown functions and the kernel functions are indicated. It is interesting 
to recognize that in (f6'2~J) in the equations of the complex unknown functions, the 
only parameter which appears in G(x) and in the integration constants, is r — k, 
which is nothing but the "length" of that part of the TBA diagram the contribution 
of which these complex functions resum. 

The energy and momentum of the model © can be easily expressed by our 
auxiliary functions due to the fact that apart from some trivial normalization factors: 

e *f V+p) ~ f fe (fe) (x + i), and e l § ^ ~ T^-xq - i). 

After some straightforward calculations ^01 EHj one gets the following result for the 
ground state energy in the continuum limit 

oo 

Eq(1) = E bulk -j J dx cosh ^"' 7 ^ log U (x) (69) 



— oo 



~ T / dx COsh ( 7r ^^ 7 ^ > ) ^ogU (x) 

— oo 

where E^uik is an overall divergent factor in the continuum and it has the form on 
the lattice: 



N 2 fe ^V 1 , rsinh[A(i(l + 2j)-xo)] 



Ebuik — -rTr - lo 

j=0 



sinh [A(i(l + 2j) + x )} 



when k is even, (70) 
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bulk — "7TT 



N 2 



(fc-3)/2 

X(2a; ) + ^ - log 



sinh [A(2z(l + j) - x ] 
sinh [A(2z(l + j) + x ] 



when is odd, 



(71) 



where 



oo 




cfg sin(gx) sinh[(r — l)q] 



(72) 



2tt q 2 cosh(g) sinh(rg) 



— oo 



From these formulas one can easily see that Ebuik is independent of the Bethe roots, 
thus independent of the dynamics of the model and so E buU , appears in the energy 
expression as a bulk energy term same for all states of the model. One can easily see 
that Eb u ik is divergent in the continuum, but after subtracting this divergent term 
from the energy the finite integral terms in ()69|) describe the finite size dependence 
of the ground state energy of the model. We note that eqs. ()62|) with the energy 
expression (}69|) are nothing else but the equations proposed by Dunning [22.J for 
the MA^ r models. However, here these equations are derived from the appropriate 
lattice regularization of these models. 

At this point we have to make an important remark, about the value of 7. 
According to Dunning' s conjecture the value of the parameter 7 can tend to 1, just 
as in the case of the sine-Gordon NLIE. Nevertheless from our derivation one can 
see, that because of the nonzero deviations of the &;-strings, 7 cannot be larger than 
1 — e, where e is the largest deviation of the fc-strings from the line lm(x) — k— 1, and 
because it is well known that e is always finite, and not infinitesimal [32!, 7 cannot 
tend to 1, and it must be kept finite. This statement can be checked numerically 
as well. Namely if one lowers the value of 7 in the numerical calculations, below 
a critical value of this parameter the numerical value of the ground state energy 
changes. In our equations we assumed that these string deviations are always less 
than 1/2 at any values of iV (or I in the continuum), that is why we restricted the 
allowed values of 7 into the (0, 1/2) interval. 
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5 Nonlinear integral equations II. 



In this section we derive a new set of nonlinear integral equations to describe the 
ground state energy of our model in finite volume, so that at the end we can glue 
it with the previous set of equations into a single set of two-component nonlin- 
ear integral equations, as it was done earlier in the S — > oo limit of the isotropic 
spin-S vertex model j2H|- Therefore let us formulate our RSOS(fc, q) model as an 
RSOS(A;', q') model with k' = r — 2 — k, q' = r — 2 — q. In this case the T-system 
(llfjjl is of the form: 

Tf \x + z)ff \x - i) = ff'\x) + f^l(x)f£i(x), (73) 

where 

ff\x) = fP(x + i(p + 1)) fit\x - i(p + 1)), (74) 

and 

fc'-l 

T {k l\x) = 0, fP{x) = J] <j>{x + i{k' - 1 - 2j)). (75) 

3=0 

The solutions of these equations are of the form 

p+i 

i=i 



where 



\ { f\x)=f^\x + i(2l-p-2)) 



e 



M P +2-2i) q(V)( x + + X )) g(fc')( x _ i(p + i)) 



Q( fc ')(x + i(2Z - p - 1)) Q( fe ')( x + z(2Z - p - 3)) ' 

(77) 



We define the following auxiliary functions: 

,f)(,) = ^^M, J = l,...,^ + l, (7E 



F/ fc,) (x) = 1 + yf >(*) = ' J v " ' V ~ ~ , J = 1, • • • , k' + 1, (79) 



, TP(x + i)fP (x-i) 

3 {X) f™{x) 



\{ k,+i \x +i)+---+Ag ; + 1 i) (x + o 



h{x) = ^ ^ f ' fc ' +1 v J -, H(x) = l + h(x), (80) 
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From (j73H77j) it follows that the auxiliary functions (j78H81|) are periodic functions 
with period 2ir and satisfy the following functional relations: 

i$2ii* + = e- (fe ' +1) + W + 2)) q^%~^ «(*), (82) 

f^U- - = ^ fP{ X - i(k> + 2)) W(x), (83) 

(X) " if^ + i(A/ + 2)) Q^(x - ik> -r) lk ' {Xh {M) 

- ~ 7777777 r; -7- T,\ 'far), (85) 

T (x-i(k + 2)) QW(x + tk + t) k 

H(x)H(x) = Y^ 1 (x), (86) 



yf'\x + i)yf'\x -i) = Y^(x)yW(x) j = 1, . . . , k' - 1, (87) 

+ + i) = Yljf\{x)H{x)H{x), (88) 



f fe ( , fc ' } (x + z)f fc ( , fc,) (a: - i) = if ' } (x + i{k' + 1)) if ' } (x - i(k! + l))^ ^). (89) 

In order to be able to derive integral equations from these functional relations one 
needs to know the positions of the zeroes and the poles of the auxiliary functions 
JZEHEU). Due to the relations flSEfl we only need the zeroes of Q {k '\x) and ff\x). 
The zeroes of Q^ k '\x) are the solutions of the Bethe Ansatz equations (j24H26j) with 
the k — > k' change. Therefore the zeroes of Q^ k '\x) form N/2 pieces of fc'-strings. 
As it was mentioned in the previous section there are deviations from the prediction 
of string hypothesis, but in the imaginary direction these deviations are always less 
than 1/2. 
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As far as the zeroes of the Tj k "* (x) transfer matrices are concerned, all Tj k ■* (x) 
can have y-fold degenerate zeroes in the "main" strip at places ±xo, ±Xo ± i, but 
these zeroes cancel from the Y-system elements (!78H79j) except the y^ ^ (x) case which 
has y-fold degenerate zeroes in the "main" strip at ±Xq. These zeroes will give the 
standard TBA source term in our final equations [H3E3|- Let us introduce 

mf'\x) = Q {k '\x - ik') (90) 

¥ k '\x) = Q {k '\x + ik') = Vf\x - 2i(r - k')) (91) 

Now we can list the strips, where the auxiliary functions are analytic and non zero 
(ANZ): 

¥ k '\x) ANZ -2(r - k') + 1/2 < Im x < 1/2, mod 2r 

^ { 2 k '\x) ANZ -1/2 < Im x < 2(r - k') + 1/2, mod 2r 

h(x),H(x) ANZ -3/2 < Im i < 0, mod 2r 

h(x),H(x) ANZ < Im x < 3/2, mod 2r (92) 

yf'\x) ANZ 0<|lmx|<l, mod 2r j = 1, . . . , k - 1, 

Yj k '\x) ANZ 0<|lmx|<e, mod 2r j = j, . . . , k e > 0, 

T^+ifc) ANZ < |Im x\ < 1 mod 2r. 

We introduce new variables by shifting the arguments of h(x),H(x) and h(x), H(x) 
by ±ry': 

a(x) = h(x — ij'), U(x) = H(x — ij') = 1 + a(x), (93) 

a(x) = H(x + irf), U(x) = H(x + i'y r ) = l + a(x), (94) 

where < 7' < 1/2 is an arbitrary real and fixed parameter. This shift is necessary 
because the original functions Tt(x), H(x) have zeroes and poles close to the real 
axis on the upper and lower half plane respectively. Due to the ANZ property of 
2fc' fc + i( x ) an d the fact that lim| a .|_ >00 log T^,^ (x) = the following relation 
holds 

00 00 
0= J dx^logf^x-i)^^- J dx^logf^ix + i)^^. (95) 



19 



From (|82j) and (J83|) one can express T^Xfa ± i) with the auxiliary functions and by 
substituting these expressions into one gets in Fourier space 

( p m (k') ( ) = e {r - k ' )q dnf{g) + e^- k ' +l -^ d>lU(q) - e^ k ~ l+ ^ SlU(q) 

4 cosh(g) 8mh[(r-k')q] 
= e^- k '^ cPlf(g) + e^- k ' +1+ ^ dW(q) - e -(^'-i-Vk ^lU{q) 
2 W ' 4 cosh(g) sinh[(r - k')q\ 

(98) 

where 

cH/(g) = d 2 / (7f >) (g) - d 2 Z {if >) (g). (99) 

After a similar procedure that has been done in the previous section one gets the 
following nonlinear integral equations in the continuum 

\ogyf\x) = (if* log Yf'^x), 

\ogyf\x) = (K*logYW) + (K*]ogYW)(x), j = 2, . . . , k' - 1 

logytfV) = + (K * \ogY^\)(x) + (K +J ' * logU)(x) 

+ (K~^' *\ogU)(x) -5n, (100) 

loga(x) = (G*\ogU)(x)-(G +2{1 -^*\ogU)(x) + (K^' *\ogY^' ] )(x) + , 

loga(x) = (G*\ogU)(x)-(G- 2 ^'U\ogU)(x) + (K^\logYP)(x) 

T rC 1 

where the "source" function on the lattice reads as 

V T N BA {x) = log ((-l) 5 tanh^ 2 (^(x - z )) tanh^ 2 (^(x + x Q ))) (101) 
which becomes the usual TBA source term in the continuum 

lim V T N BA (x) = -ml cosh (^-x) . (102) 

N—*oo V 2 / 

The kernel function K(x) and S are the same as in the previous section and 

oo oo 

e iqx ginh^ _ _ 2)g] f dq e iqx sinh(A;g) 



G(x) = / -± 



27r 2 cosh(g) sinh[(r — k' — l)q\ J 2n 2 cosh(g) sinh[(/c + l)g] ' 

— oo 

(103) 
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The correct choice of the value of 5 is zero for the ground state of the A4A^ r . 
In the TBA language the complex auxiliary functions a{x) and a(x) resum the 
contributions of those Y-system elements whose index is larger than k'. (See figure 

E3>0 

Due to the identity (x) = Y^,(x), which follows from (J3Hj) the energy (0) 
of the continuum model can be expressed by a gauge invariant Y-system element of 
the RSOS(fc',g') model: 

oo 

E Q (l) = E bulk -j J cfecosh (H) logics). (104) 

— oo 

Equations (J52j) and (|100j) are different descriptions of the ground state energy of the 
massive 4>^ ac |j perturbation of the SU(2)k X SU{2)k> / 'SU(2)k+k' model. In this 
paper one of our main purposes is to describe the finite size dependence of the ground 
state energy of this model with a two-component nonlinear integral equations. As 
we will see in the next section following the idea of j22j we are able to construct 
such two component nonlinear integral equations by combining (J52j) and (jlOOj) . 

6 The two component nonlinear integral equa- 
tions 

In this section we construct two component nonlinear integral equations from (pT2*|) 
and ([100)1 . The complex auxiliary functions ()44M5|) and (jSPHS!)) are gauge dependent 
and have a finite continuum limit only in an appropriate gauge, but the Y-system 
elements (j29H31j) are gauge invariant and have finite continuum limit. One can see 
from (|50|) and ()86|) that some gauge invariant Y-system elements can be expressed 
by these non-gauge invariant auxiliary functions: 

U (x + i^Uoix - i 7 ) = Y^ k) (x), (105) 
U(x + ij')U(x - ii) = Y^\(x). (106) 
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In our case because of the equivalence of the Y-systems of an RSOS(A;,g) and 
KSOS(k',q') model (|35|) the following relation will also be true 

U {x + ii)U {x - z 7 ) = Y$'\x). (107) 



U{x + ii)U{x - ii) = Y^\{x). (108) 

Substituting ()108j) into the equations for ao(x),ao(x) in (|H2^I and ()107j) into the 
equations for a(x),d(x) in (jlOOj) and deforming the integration contour appropri- 
ately, we get the two component nonlinear integral equations for the ground state 
of our model (5 = 0). The equations are as follows: 



loga (x) = V{x-ij) + (G*logU )(x)-(G +2 ^ *logU )(x) 
+ (K + ^'-^ * log U) (x) + (if-^+T) * log U) (x) + % 



71 



r — k' 

logao(x) = P(a: + « 7 ) + (G*logf/ )(a;)-(G'- 2(1 - 7) *\ogU )(x) 

+ (if+fr'+r) * logU)(x) + (iT +(7 " 7,) * logC/)(x) - i-^-r, 

r — k 

\oga(x) = {G*logU){x)-{G +2{1 - y) *logU){x) (109) 

+ # l O gf/ )(x) + (K-^' + ^ * \ O gU )(x)+i-^-, 

k + 1 

loga(x) = (G*log[7)(x) - (dr 2 ^') *\ogU)(x) 

+ (K+^+t') * log I/o) (x) + (K +{7 '~ 7) * log C/ ) (x) - i 
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k + r 

U (x) = l + a (x), Uo(x) = 1 + a (x), U(x) — l + a(x), U(x) = 1 + a{x). 

where < 7 < 1/2, < 7' < 1/2 are arbitrary fixed real parameters. Equations 
(I109J) are graphically represented in figure Eh- 



O — O : — • — O O — O a 

1 2 k-iy \k k+l r-4 r-3/ 



a(x), a(x), G(x) a (x),a (x),G(x) a{x),a{x),G(x) a (x) , a (x) , G(x) 

Figure 3: Graphical representation of the two component nonlinear integral 
equations of the o.) MA^ r and b.) FSSG model 

The energy expression is the same as (jEHJ)- This is a closed set of equations for 
four complex unknown functions "a(x), a(x), a (x), a (x)" , but a(x) is the complex 
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conjugate of a{x) and a (x) is the complex conjugate of a (x), therefore only four 
real unknown functions describe the model. 

For the special case of r = 2k + 2, from (jlUUJI one can derive a different set 
of nonlinear integral equations to describe the finite size dependence of the ground 
state energy of our model. In this case k = k! and thus Yy\{x) = Y^_ x (x) = 
U{x + i^')U(x — 27'). Replacing this into the equation of y^ \x) in (jlOOj) 

and deforming the integration contour appropriately we get the following equa- 
tions 

\ogy (x) = V(x) + 2(K + " / ' * log 17) (2;) + 2{K^' *\ogU){x), 
loga(x) = (G*logU)(x) - (G +2(1 - 7,) *logU)(x) 

+ {K-i 1 *i og y )(x) + ^ I , (110) 

loga(x) = (G*\ogU){x) - (CT 2(W) *\ogU){x) 



+ (K+i' *logY )(x) 



k+V 

U(x) = l + a(x), U(x) = 1 + a(x), 
where we introduced the notations: 

y (x)=y ( *\x), Y Q {x)=Yt k \x). (Ill) 

Equations ()110j) are graphically represented in figure EJ 



3 O O 

K k+l 2k+2 



a(x),a(x),G(x) a(x), a(x), G(x) 

Figure 4: Graphical representation of the NLIEs of the MA^ 2k+2 model. 

In this case the ground state energy of the model is of the form 

00 

Eq(1) = E bulk -— da; cosh J \ogY (x). (112) 

—00 

Equations ()110|) contain only three real unknown functions because yo(x) and Y (x) 
are real and a(x) is the complex conjugate of a(x). 
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It can be easily seen that the k — > oo limit of equations (|110|) are the same as 
the ones which were proposed earlier j2S| for the 0(4) NLS model. For finite k 
these equations describe the M-A^ 2k+2 model, the k — > oo limit of which is the 0(4) 
nonlinear cr-model fTf]. 



7 The test of the equations 

In this section we will make some analytical and numerical tests on our equations 
First we calculate the UV central charge of the model using the equations 
(jl()9j) . Using the standard method of refs. [2HEEJ the energy in the conformal limit 
(/ — > 0) can be expressed by the dilogarithm functions. The energy in the conformal 
limit is of the form: 

E o(l) - { £+( a +(°o)) + L + (a + (oo)) + L + (a +(oo)) + L+(a 0+ (oo)) 



irl x 

-L + (a + (-oo)) - L+(a+(-oo)) - L+(a +(-oo)) - L + (a +(-oo))^ (113) 
+i 2(7^k) ^ " ^o+(+oo)) - (L4 0+ (-oo) - L4 0+ (-oo)) } 

^ ( M+(+oo) ~ L4 +(+°°)) - - £*+(-°°)) } 

where 

lA+(x) = log(l + a+(x)), L4 + (x) = log(l + a+(x)), (114) 
lA 0+ (x) = log(l + a +{x)) l~A 0+ (x) =log(l + a 0+ (x)), 

L + (z) is defined by the integral 

i+W = i/ d J!2i(l±f)_ T ^y L+ (x) = l(^-), (115) 



2 J \ X 1 + X J \1 + X 



and L(x) is Roger's dilogarithm function. The functions a+(x), a+(x), a 0+ (x), a 0+ (x) 
denote the kink functions corresponding to a(x), a(x), a (x), a (x) respectively. The 
limits of these kink functions at infinity are as follows: 

sin ( (fc, + 1)7r 

a + (oo) = (a+(oo))* = e*^, a+(-oo) = (a+(-oo))* = e 4 ^ 2 ^ . 

(116) 
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(fc+i) sin ( 

a 0+ (oo) = a 0+ (oo) = 0, a 0+ (-oo) = (a 0+ (-oo))* = e*^~ , *\ , (117) 

sm (2) 

where the * denotes the complex conjugation. From these using some dilogarithm 
identities jHE], one gets the ground state energy in the I — > limit 



O I 

Where qyy(/c,r) is the effective UV central charge of the model, the value of which 
according to our equations is 

3k ( 2 (k + 2) \ 

Cuv ^ r) = —2\ l ~ (119) 

which agrees with the UV central charge of the A4A^ r model. So we have checked 
our equations (jl()9j) analytically in the small I limit in leading order. Next we will 
check our equations in the large I limit analytically in leading order. The equations 
(\109\i can be solved iteratively in the large / regime. After some easy calculations 
one gets that the ground state energy is of the form 

E (l) = E {1) +0(e- 2ml ), (120) 

where 

oo 

Em = - 4 <=» (7^1) » s (kh) S de r, cosm e ' m "" m ' (121) 

— oo 

which agrees perfectly with TBA results obtained for M.A^ r models 0. Note, that 
this expression is in accordance with the expectation for the low-temperature be- 
haviour of the kink system with the vacuum structure described in the introduction. 
In (J121)) the factor in front of the integral expression takes into account the statistics 
of interkink "colouring" with the ® ^4fc'+i adjacency structure. Therefore this 
limit is in agreement with the (RSOS)^ x (RSOS)fc/ scattering theory of ref.0. 

So far we have made some analytical tests on our equations ()109|) in the small / 
and large / regimes. We also made numerical calculations to check the correctness 
of our equations (jl()9j) . Therefore for different values of the parameters k,r and /, 



25 



we solved numerically our eqs. (jl09J) . the corresponding TBA equations [H], and 
Dunning's "half TBA-NLIE" type of equations [22] • In every case all the three 
type of equations gave the same numerical results, which convinced us about the 
correctness of our equations (jl09|) . 

8 Two component nonlinear integral equations for 
the fractional supersymmetric sine-Gordon mod- 
els 

In this section we propose two component nonlinear integral equations for the frac- 
tional supersymmetric sine-Gordon models. One can recognize that taking the 
r — > oo limit of eqs. (105), one obtains two component nonlinear integral equa- 
tions for the current- current perturbation of the SU(2)k WZW model, which is a 
fractional supersymmetric sine-Gordon model at a special value of the coupling. The 
equations take the form 

log a (x) = V(x - ij) + (F * log U ) (x) - (F +2(1 " 7) * log U ) (x) 

+ (K+^-t) * log[/)(» + (K-^' + ^ * log [7)0), 

logao(x) = V(x + ij) + (F*\ogUo)(x)-(F- 2 ^*\ogUo)(x) 

+ (K+« + ^ * log U) (x) + (K^-^ * log U) (x) , 

loga(z) = (G*\ogU){x) - (G +2(W) * log 17) (a;) (122) 

+ * \ogU )(x) + (K-W-^ * \ogU ){x)+t-^-, 

K + 1 

loga(x) = (G* log 17) (a) - (G-^ 1 -^ *\ogU){x) 



+ (K + ^' ] * log U ) {x) + (K + ^-^ * log U ) {x)-i- 



k_+ V 

U Q (x) = 1 + a (x), Uq(x) = 1 + a (x), U(x) = l + a(x), U(x) = 1 + a(x). 

where where < 7 < 1/2, < 7' < 1/2 are arbitrary fixed real parameters, G(x) is 
the same as in (100), but 

00 

f dq e -\ q \- iqx . , 

Fix) = / — r-T, 123 

K J J 2tt 2 cosh g ' V ; 
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is the r — ► oo limit of G(x), and the energy formula agrees with (|69|). One can see 
in eqs. (1122jl . that the variables a(x), a(x) resum the contributions of those nodes of 
the corresponding Dynkin-diagram, whose index are less than k. The lattice version 
of these equations correspond to the case of the isotropic spin-fc/2 vertex model. 
So far, all of the two component nonlinear integral equations we produced, were 
the results of derivations. Now we will propose two component nonlinear integral 
equations for the fractional supersymmetric sine-Gordon models, the lattice version 
of which will correspond to the case of the anisotropic higher spin vertex model. 

One can see in eqs. (|122|) that those components of the equations, which have 
log a(x) and \oga{x) on the left hand side, are independent of r (the entire length 
of the Dynkin-diagram) and the only parameter which appears in them is the k, the 
length of the resummed part of the Dynkin diagram. From this one can conclude, 
that the form of these two components of the equations remain the same, if the right 
hand side of the Dynkin-diagram is changed from A r _ s to D r . Therefore deforming 
the kernel function F(x) of eqs. ()122|) to G(x), as it was done in [22J, we obtain 
two component nonlinear integral equations for the FSSG models. The proposed 
equations are as follows 

logao(x) = V (x-i-f) + (G* log U )(x)-(G +2{1 -^ *\ogU )(x) 

+ (K +i " f '-" f) * \ogU)(x) + (tf-W+r) * \ogU){x), 

loga (x) = V(x + i-f) + (G*\ogU )(x)-(G- 2 ^ ) *\ O gU )(x) 

+ (K+W+rt * log U) (x) + (if * log U) (x) , 

loga(x) = (G*logU)(x) - (G +2(1 - 7 ' } *logU){x) (124) 

+ (K+b-rt * log U ) (x) + (K-(t'+t) + i og f7 ) (x)+i. 
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k + r 

hga(x) = (G*\ogU)(x) - (G~ 2(W) *\ogU)(x) 



+ (K + ^+^ * log U ) (x) + (K+(t'-t) * log U ) (x)-i-. 
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k + r 

U (x) = l + a (x), Uq(x) = 1 + ao(x), U(x) = l + a(x), U(x) — 1 + a(x). 



where kernel functions are the same as in and the energy expression is also 

the same (joTJj) . These equations describe the ground state energy of the fractional 
supersymmetric sine-Gordon models, where the relations between the f3 2 coupling 
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of the model, and the parameters r and k of the eqs. (j!24|) are 

(3 2 r-k 



(125) 



8n r k 

For the graphical representation of the equations ()124j) see figure From the 
eqs. fll24j) after some usual straightforward calculations one can easily obtain the 

3fc 

k+2' 

behavior of the energy one obtains 



expected value of the central charge cuy = -Mz, and for the leading order infrared 



oo 

E (l) ~ -4 cos ( -^-) ^ [ d6 cosh(#) e - mZcosh(e) , (126) 
\k + 2/ 27r y 

— oo 

which agrees perfectly with TBA results. 

We mention, that eqs. (|124|) could be obtained by the following formal train 
of thought: the FSSG models were first proposed as the models obtained by "un- 
restricting" the S-matrices describing M.A^ r models ilj. Eqs. ()109|) describe the 
finite size effects in this model, and so we can formally "release" the restriction by 
removing the r dependent imaginary constant terms from (jlOfJJ) . and thus we obtain 
eqs. ([124)1 . Then by formally "releasing" the remaining restriction in eqs. ()124|) 
by removing the remaining /c-dependent imaginary constant, one obtains the earlier 
proposed two component nonlinear integral equations of the SS-model |2Hj . 

We also made numerical calculations to check the correctness of our conjectured 
equations ()124j) . Therefore for different values of the parameters k,r and I, we 
solved numerically our eqs. ()124)) . the corresponding TBA equations [J], and the 
Dunning's "half TBA-NLIE" type of equations [22]- in every case all the three 
type of equations gave the same numerical results, which convinced us about the 
correctness of our equations (J124)) . 



9 The equations for finite solvable lattices 

In this section we briefly summarize how the two component nonlinear integral 
equations of the previous sections have to be modified for the case of solvable lattice 
models. In the context of lattice models our purpose is to calculate exactly the 
eigenvalues of the transfer matrices (fT4"j) (with xq = 0) for finite N, because by 
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taking the first derivative of these transfer matrices at certain points one can define 
integrable spin Hamiltonians from them ITTI] . The derivation described in the 
previous sections can be achieved for the largest eigenvalue of the transfer matrix of 
the critical RSOS(/c, q) model. The equations are the same as (jl(J9|) with the change 
of the source function T>(x) to 



T>n(x) = iN arctan 



sinh 



TT[X + i) 



2 

The largest eigenvalue of the transfer matrix is given by 



(127) 



oo 

log (a;, iV) = log A$ ulk (x, N) + \ [ dx> lo f^ X ') (128) 

4 _4 COsh M»-^+*7? j 



+ - [ dx' 1Ogf/0 ^ 



4 J cosh ( <x-x'-n) 
— oo y £ 

where log A^l ulk (x, N) is the term proportional to N. By taking the r — > oo limit 
of the equations of the critical RSOS(fc,g) lattice model, we obtain eqs. ()122j) with 
the T>(x) — > T>n(x) change of the source function, which describe the A" dependence 
of largest eigenvalue of the transfer matrix T k (x) of the isotropical spin-fc/2 vertex 
model. The equations ()124|) with the V(x) — > T>n(x) change are though to describe 
the the finite size effects in the anisotropic higher spin vertex model at the spin value 
S = k/2. The 1/N corrections to the largest eigenvalue of the transfer matrices can 
be analytically calculated and in every case the well-known results were obtained. 



10 Summary and Conclusions 

From the inhomogeneous RSOS(A;,g) lattice model considered on the critical line 
separating regime III and IV, we derived two different sets of nonlinear integral 
equations describing the finite-size effects in the lattice model. The appropri- 
ate continuum limit of the equations describe the finite-size dependence of the 
ground state energy of the massive 0^ ac y perturbation of the coset SU(2)k x 
SU(2) r ^k-2/ SU(2) r -2- These two sets of nonlinear integral equations are based 
on the two different formulation of the same model. One set of equations is de- 
rived in the RSOS(fc,g) model, and the other set of equations is derived in the 



29 



RSOS(r — k — 2, r — q — 2) model, using their equivalence. These two formulations 
correspond to two different representation of the solution of the same (apart from a 
trivial gauge transformation) T-system. These two different representations of the 
solutions are characterized by the solutions of two different Bethe Ansatz equations. 

Then we glued these two sets of nonlinear integral equations together and ob- 
tained two component nonlinear integral equations for the finite size effects. Taking 
the r — > oo limit of our equations we obtained two component nonlinear integral 
equations corresponding to the case of the isotropic six- vertex model, and the contin- 
uum limit of which describe the finite-size effects in the current- current perturbation 
of the SU(2)k WZW model. Modifying the kernel functions of these equations we 
proposed two component nonlinear integral equations for the fractional supersym- 
metric sine-Gordon models. The lattice version of these equations are thought to 
describe the finite-size effects in the anisotropic higher spin vertex model. 

We performed analytical and numerical tests of our equations and all the results 
of these tests made us confident that our equations are correct. The advantage 
of our equations is that the number of unknown functions are minimal, and in 
the case of the FSSG models the equations can be defined for all real values of 
the coupling constant. Another advance of having such two component nonlinear 
integral equations appears in the investigation of the excited state problem of the 
models, where only two components have to be treated, instead of the large number 
of components of the TBA equations. The treatment of excited states in our models 
is left for future investigation. 
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